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ABSTRACT 

We analyze the surface structure of the planet host star CoRoT-2a using a consistent model for both the 'global' (i.e., 
rotationally modulated) lightcurve and the transit lightcurves, using data provided by the CoRoT mission. Selecting 
a time interval covering two stellar rotations and six transits of the planetary companion CoRoT-2b, we adopt a 
'strip' model of the surface to reproduce the photometric modulation inside and outside the transits simultaneously. 
Our reconstructions show that it is possible to achieve appropriate fits for the entire sub-interval using a low-resolution 
surface model with 36 strips. The surface reconstructions indicate that the brightness on the eclipsed section of the stellar 
surface is (6 ± 1) % lower than the average brightness of the remaining surface. This result suggests a concentration of 
stellar activity in a band around the stellar equator similar to the behavior observed on the Sun. 

Key words, planetary systems - techniques: photometry - stars: activity - stars: starspots - stars: individual: CoRoT-2a 
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1. Introduction 

Astronomers have long been interested in the surface struc- 
ture of active stars and their evolution. Yet, the surfaces of 
stars other than the Sun can hardly be resolved directly, 
so that indirect techniques must be used to obtain an im- 
age of the surface. One such technique is Doppler imag- 
ing (Vogt & Penrod 1983), which requires a dense series 
of high resolution spectra and stellar rotation velocities of 
vsm{i) > 20 km/s (compared to Vgq ~ 2 km/s for the 
Sun). Alternatively, lightcurves also yield information on 
stellar surface structures and can usually be obtained at low 
observational cost. However, photometry provides less in- 
formation and the problem of lightcurve inversion is known 
to be notoriously ill-posed. 

Since the launch of CoRoT in 2006, an increasing 
amount of high quaUty space-based photometry has be- 
come available. Without the limitations the atmosphere 
and the day-night cycle impose on ground based observa- 
tories, CoRoT is able to provide photometry with unprece- 
dented temporal coverage and cadence, which is enormously 
interesting in the context of stellar activity and surface re- 
construction. 

In the course of the CoRoT planet hunting project, the 
giant planet CoRoT-2b (Alonso et al. 2008) was detected. 
The host star of this planet, CoRoT-2a, is solar like in mass 
and radius, but rotates approximately four times faster 
than the Sun and is considerably more active. The planet 
orbits its host star approximately three times per stellar 
rotation and, during its passage across the stellar disk, acts 
as a shutter scanning the surface of the star along a well 
defined latitudinal band. As the 'local' surface structure is 
imprinted on the transit profiles (Wolter et al. 2009; Czesla 



et al. 2009), they can be used to partially resolve the am- 
biguity of the lightcurve inversion problem. 

While Lanza et al. (2009) used the 'global' lightcurve 
of the host star to reconstruct its surface inhomogeneities, 
without considering the transits, Wolter et al. (2009) con- 
centrated on a single transit lightcurve to reconstruct a frac- 
tion of the surface, neglecting the 'global' lightcurve. In this 
work, we combine and refine these approaches to present 
a reconstruction which simultaneously describes both the 
overall lightcurve and the transits during two stellar rota- 
tions. 



2. Observations and data reduction 

Alonso et al. (2008) discovered the planet CoRoT-2b using 
photometric data provided by the CoRoT mission (for a 
detailed description, see Auvergne et al. 2009). The planet 
was detected in the field observed diiring the first long run 
carried out between May 16"^ and Oct. IS^'' 2007. The de- 
fault sampling rate of CoRoT photometry is 1/512 s^^. 
The CoRoT-2 lightcurve was observed at this rate only for 
the first five days, after which the transits were detected 
and the satellite switched to alarm-mode, now taking data 
every 32 s. The light collected by the CoRoT telescope 
is dispersed using a prism and recorded by a CCD chip. 
Individual sources are separated by a photometric mask, 
which also defines three broad-band channels (nominally 
red, green, and blue). Currently, there is no appropriate 
calibration available for these channels, so that it is unfea- 
sible to use the color information in this work. The signal 
obtained by summing up the individual channels, often re- 
ferred to as 'white light', corresponds to an optical measure- 
ment with a filter transmission maximum in the red wave- 
length region (Auvergne et al. 2009). Accordingly, Lanza 
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Table 1. Stellar/planetary parameters of CoRoT-2a/b. 



Star ° 


Value zb Error 


Ref.' 


S 


(4.522 ± 0.024) d 


L09 


s 


4.57 d 




spectral type 


G7V 


BOS 


PlaTic.+ 

r^ianet 


Value i Error 


iiei. 


Pp 


(1.7429964 ± 0.0000017) d 


AOS 


n [BjD] 


(2 454 237.53362 ± 0.00014) d 


AOS 


i 


(87.7 ± 0.2)° 


C09 


Rp/Rs 


(0.172 ± 0.001) 


C09 


a/Rs 


(6.70 ± 0.03) 


AOS 


Ua, Ub 


(0.41 ±0.03), (0.06 ±0.03) 


AOS 



" Ps - stellar rotation period, Pg - stellar rotation period used 
for the observation interval analyzed in this paper (see Sect. 3.4). 

' taken from Lanza et al. (2009) [L09], Alonso et al. (200S) 
[AOS], Bouchy et al. (200S) [BOS], or Czesla et al. (2009) [C09] 

° Pp - orbital period, Tc - central time of first transit, i - or- 
bital inclination, Rp,Rs - planetary and stellar radii, a - semi 
major axis of planetary orbit, Ua,ui, - linear and quadratic limb 
darkening coefficients. 



et al. (2009) assume an isophotal wavelength of 700 nm 
for their pass-band. The CoRoT data undergo a standard 
pipeline processing, during which data points that are sig- 
nificantly affected by known events, as for example the 
South Atlantic Anomaly (SAA), are flagged, so that they 
can be removed from the lightcurve. 

The host star CoRoT-2 has a spectral type of G7V with 
an optical companion at a distance of approximately 4.3" 
(2MASS, Skrutskie et al. 2006), too close to be resolved by 
CoRoT. According to Alonso et al. (2008) the secondary 
contributes a constant fraction of (5.6 ± 0.3) % to the 
total CoRoT-measured flux. In Table 1 we list the sys- 
tem parameters of CoRoT-2a/b, which are used throughout 
our analysis. CoRoT-2b's orbital period of « 1.74 days is 
about a third of CoRoT-2a's rotation period. Hence, the 
almost continuous CoRoT data sample of 142 days covers 
about 30 stellar rotations and more than 80 transits. The 
lightcurve shows signatures of strong stellar activity and 
substantial rotational modulation (Lanza et al. 2009). We 
use the same CoRoT raw data reduction procedures as de- 
scribed in Czesla et al. (2009, Sect. 2). 



3. Analysis 

3.1. Modeling approach 

The measurements of the Rossiter-McLaughlin effect by 
Bouchy et al. (2008) suggest that the rotation axis of 
the host star and the planet's orbit normal are approxi- 
mately co-aligned (A = 7.4 ± 4.5°). A represents the mis- 
alignment angle projected on the plane of the sky, and its 
value strongly favors aligned orbital and rotational axes, 
even though it does not prove it. Further support for a 
co-aligned geometry comes from the following argument: 
Comparing the measured wsin(j) = 11.85 km/s with a cal- 
culated equatorial velocity of v^q = 2-KRg/Pg w 10 km/s 
derived with the theoretically obtained value Rg — 0.9 • Rq 
(Alonso et al. 2008) also favors sin(i) « 1. 



Therefore, the planet always eclipses the same low- 
latitude band between 6 and 26 degrees. The transits sep- 
arate the stellar surface into two observationally distinct 
regions, i.e., a region eclipsed by CoRoT-2b and another 
region that is not. In the case of CoRoT-2a the eclipsed 
section covers ~ 21 % of the stellar disk corresponding to 
« 17.3 % of its surface. The time-resolved planet migra- 
tion across the visible stellar disk sequentially covers and 
uncovers surface fractions, so that the brightness profile of 
the underlying stellar surface is imprinted on the transit 
lightcurve. 

For our modeling, we separate the surface into the 
eclipsed and the non-eclipsed section, which are both fur- 
ther subdivided into equally sized, longitudinal bins or 
'strips' as demonstrated in Fig. 1. Let iVg be the number 
of bins in the eclipsed section and Nn be the number of 
non-eclipsed strips. As is apparent from Fig. 1, and 7V„ 
need not be the same. Altogether, we have Ntot = Ne + Nn 
bins enumerated by some index j. To each of these surface 
bins a brightness bj is assigned, with which it contributes 
to the total (surface) fiux of the star. Now let Vji denote 
the visibility of the j-th bin at time U. The visibility is 
modified in response to both a change in the viewing ge- 
ometry caused by the stellar rotation and a transit of the 
planet. The modeled flux fmod,i at time U is then given by 
the expression 

Nt^t 

fmod,i — ^ ^jibj ■ (1) 

We determine the unknown brightnesses, bj, by comparing 
fmod,i to a set of Mc CoRoT flux measurements using a 
specifically weighted version of the x^-statistics: 

Xm = 2^ 2 ■ (2) 

i=l 

Xm differs from by a weighting factor, Wj, which we 
choose to be 10 for lightcurve points in transits and 1 oth- 
erwise; in this way the global lightcurve and the transits 
are given about the same priority in the minimization pro- 
cess. Error bars for the individual photometric measure- 
ments were estimated from the data point distribution in 
the lightcurve, and the same value of tr = 1000 e^/32 s 
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Fig. 1. Our model geometry using 12 longitudinal strips for the 
non-eclipsed and 24 strips for the eclipsed section, respectively. 



K. F. Huber et al.: A planetary eclipse map of CoRoT-2a 



3 



(= 1.4 X 10~^ after lightcurve normalization) was used for 
all points. 

In our modeling we currently exclude surface structures 
with a limb-angle dependent contrast; this particularly 
refers to solar-like faculae, for which Lanza et al. (2009) 
find no evidence in their analysis. The planet, CoRoT-2b, is 
modeled as a dark sphere without any thermal or reflected 
emission. This approximation is justified by the findings 
of Alonso et al. (2009), who report a detection of the sec- 
ondary transit with a depth of (0.006 ± 0.002) %, which is 
negligible in our analysis. 

The actual fit is carried out using a (non-gradient) 
Nelder-Mead simplex algorithm (e.g. Press et al. 1992). In 
our model all strips are mutually independent, and as we 
define only a rather coarse strip subdivision for the non- 
eclipsed section of the surface and since the eclipsed section 
is thoroughly covered by the transits, no further regulariza- 
tion is necessary. 

3.2. Fit parameter space - parametrization, restrictions, and 

interpretation 

Our fit space has a total of Ntot = N^ + Nn dimensions, and 
the associated parameters are the brightnesses, bj=i,,N^^^. 
The most obvious choice of fit parameters are the bright- 
nesses themselves. Yet, using a slightly different definition 
in our algorithm provides some advantages. Instead of us- 
ing the brightness of the global strips in our fits, we replace 
them with a weighted sum of the brightnesses of the eclipsed 
and the non-eclipsed strips covering the same longitudes. 
This quantity z is a, measure of the total flux emitted from 
all strips enclosed within a certain longitude range and, 
therefore, represents the level of the global lightcurve inde- 
pendent of how the brightness is distributed among the in- 
dividual strips contributing to the sum. Without any transit 
observations the distribution of flux among the individual 
contributors could hardly be further restricted, because lat- 
itudinal information could not be recovered. Thus, we use 
the tuple {bi=i,,N^,Zk=i..N„) for our reconstructions, where 
z is defined by 

^ S<IJO+IJ 

Zk = bNe+k H bs . (3) 

c - q ^ 

In Eq. 3 bN^+k denotes the brightness of the k-th global 
strip, q is defined by N^/Nn (the factor by which the 
ecUpsed section is oversampled compared to the non- 
eclipsed section), the index range qn < s < qo + q enu- 
merates all eclipsed strips covering the same longitudes as 
the global strip referred to by b^^+k, and c is a scaling fac- 
tor accounting for the size difference between the eclipsed 
and the non-eclipsed section. 

The practical advantage of using z instead of the bright- 
ness values themselves lies in the parameter interdepen- 
dence. Assume a fit algorithm adjusts the structure of a 
transit lightcurve using the eclipsed strips; every modifica- 
tion of their brightness causes a modification of the global 
lightcurve level which must possibly be compensated by 
an appropriate adjustment of the global strip's brightness. 
Such an adjustment is inherent in the definition of z, so that 
bi=i..Nc and Zk=i..N„ become largely independent quanti- 
ties. In our fits we use c = 5, which roughly corresponds to 
the ratio of disk area covered by global and eclipsed strips. 



To normalize the observed CoRoT-2a lightcurve, we di- 
vided all measurements by the largest fiux value in our 
lightcurve so that < (normalized flux) < 1. The matrix 
Vji in Eq. 1 is normalized according to 

Ntat 

Vji = 1 for all i , 

which yields fmod,i = 1 for bj — 1, i.e., a constant model 
lightcurve at level 1. In a flrst, tentative interpretation a 
surface element with the brightness 1 corresponds to a pho- 
tospheric element free of any spots. Yet, this is only cor- 
rect as long as we assume that the largest observed flux in 
the lightcurve, indeed, represents the 'spot-cleaned' photo- 
spheric luminosity. As CoRoT-2a is, however, a very active 
star, it seems probable that polar spots persist on its sur- 
face. Moreover, it seems likely that lower latitude structures 
cover a fraction of the stellar surface even if the lightcurve is 
at maximum. For this reason, individual surface elements 
(strips) may be brighter than the 'average surface' dur- 
ing the maximum observed flux. While such information 
could not be recovered if no transits were observed, indi- 
vidual surface regions eclipsed by the planet can conceiv- 
ably be brighter than the 'global' photosphere seen during 
lightcurve maximum. Therefore, we do not exclude strips 
with brightness values larger than 1 in our fits, i.e., we do 
not fix the photospheric brightness; this results in bright- 
nesses greater than 1 for individual strips (e.g. Fig. 5). The 
only parameter space restriction applied during our fits is 
that the brightness must be positive. 

3.3. Which part of the lightcurve should be used? 

In order to derive a meaningful model we need to select 

a time span, which is both long enough to provide an 
appropriate coverage of the surface, and short enough to 
minimize the effects of surface evolution; the latter, while 
doubtlessly present, appears slow compared to the stellar 
rotation period. Lanza et al. (2009) give typical lifetimes of 
55 d (w 12 rotations) for active regions and 20 — 30 d for 
some individual spots. In our analysis, we use the time span 
ranging from phase 1.85 through 3.85 (BJD = 2 454 245.988 
to BJD = 2454255.128, BJD = Barycentric Julian Date), 
which covers 6 transits and shows only little variation in 
the global lightcurve. The data are re-binned using a bin- 
size of 128 s for the transit covered periods and 2 016 s for 
the remaining lightcurve. 

Our binning approach has to take into account interrup- 
tions of the lightcurve due to data drop outs (for instance 
due to the South Atlantic Anomaly) and, of c;ourse, has to 
account for the change in bin size when a transit period be- 
gins or ends. Moreover, the CoRoT-2 lightcurve is sampled 
at two different rates (1/512 s""'^ and 1/32 s""'^), which does, 
however, not impose a problem during the time span under 
consideration here. To obtain the binned curve, we average 
all flux values comprised by a bin and place the resulting 
value at the barycenter of the contributing time stamps. 
To compute the error, we divide the standard deviation for 
individual points by the square root of the number of av- 
eraged points. With this approach we (typically) obtain an 
error of 7 x 10~^ for in-transit points and 1.8 x 10~^ for 
out-of-transit points. 

In Fig. 2 we demonstrate the coverage of the eclipsed 
surface section by these 6 transits within the selected phase 
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interval. A single rotation phase including three transits 
provides only a very inhomogeneous 'scan' of the eclipsed 
surface due to limb darkening, projection geometry, and 
the distribution of transit intervals (cf. Fig. 2). As a transit 
occurs every « 0.4 stellar rotations, a homogeneous cov- 
erage of one full rotation is achieved using five transits. 
Nonetheless, we decide to use an integer number of stel- 
lar rotations and use six transits with the last one showing 
virtually the same part of the eclipsed surface as the first. 

3.4. Surface evolution, rotation period, and model limits 

Although the lightcurve of CoRoT-2a shows remarkably pe- 
riodic minima and maxima, the rotation period of the star 
is not exactly known. Using a Lomb-Scargle periodogram, 
Lanza et al. (2009) find a rotation period of (4.52 ± 0.14) d 
for the star, which is further refined in the course of their 
surface modeling. Assuming that the longitudinal migra- 
tion of the active longitudes should be minimal, Lanza et al. 
pin down the stellar rotation period to 4.5221 d. While this 
rotation period minimizes the migration of the active lon- 
gitudes, it results in individual spots showing a retrograde 
migration with an apparent angular velocity « 1.3 % lower 
than the stellar rotation. 

For our modeling we determine an 'effective' period, rep- 
resenting the rotation period of the dominating surface fea- 
tures we are mainly interested in. In our approach, we use 
the selected part of the lightcurve, remove the transits, and 
fold the remaining lightcurve back at a number of periods 
between 4.4 d and 4.7 d. The best match is obtained using 
a period of 4.57 d, which also results in the best fits of our 
models; therefore, we will use it throughout our analysis. 
This period is also in agreement with the values given by 
Lanza et al. (2009) assuming a rotation period of 4.5221 d 
and a mean retrograde migration 'slowing down' the spots 
by 1.06 %. Note, however, that changes of the rotation pe- 
riod on this scale do not result in significantly different 
surface reconstructions. 

Even though we identified a lightcurve interval with rel- 
atively weak surface evolution, and refine the rotation pe- 
riod to account for some evolutionary effects, there is still 
a remaining modulation. The presence of this modulation 
imposes a principal limit on the fit quality that can be 
achieved by adopting a static model to the lightcurve, be- 
cause both stellar rotations have to be described by the 
same model. To estimate this limit, again for the global 
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Fig. 2. Visibility of the eclipsed stellar surface during the tran- 
sits in the selected observation interval. A low visibility means 
that a stellar feature at the corresponding longitude has a low 
impact on the transit profiles. 



lightcurve only, we estimated the quantity 

i=l ^ ^ 

Here, f(j>i) is the normalized fiux in the i-th phase bin, 
f{pi + 1.0) is the fiux measured at the same phase dur- 
ing the next stellar rotation, and the sum stretches over 
all phases pertaining to the first rotation. f{pi + 1.0) was 
obtained by interpolation, because the phase sampling is 
not exactly the same in both rotations. Since Ui « cr, the 
best conceivable common model with respect to phase 
point Pi is given by {f{pi) + f{pi + 1.0))/2, and, therefore, 
the sum in Eq. 4 estimates the contributions induced by 
surface evolution for each point. If there was no surface evo- 
lution, the expression in Eq. 4 would equate to 0.5, because 
statistical errors are, of course, still present. Therefore, a 
limit of « 14.2 per (global) lightcurve point will not be 
overcome by any static model. Equivalently, the expecta- 
tion value, < A//2 >, for the fiux deviation from the best 
model equates to < A//2 >=< (/(k) - f{p, + 1.0))/2 >= 
5.6 X 10^^ and cannot be surpassed. 

3.5. Model resolution 

The parameters A^e f^nd specify the model resolution 
of the eclipsed and non-eclipsed sections. An appropriate 
choice of these parameters balances fit quality and model 
ambiguity; this way the largest possible amount of infor- 
mation can be extracted. 

In order to find the optimal value for the number of non- 
eclipsed strips, we carry out fits to only the global lightcurve 
using an increasing number of global strips. Starting with 
only 4 strips, we find the reduced value, Xjii to decrease 
rapidly until 8 strips are used. From this point, x\ responds 
only weakly to an increase of the strip number, but still 
decreases. Using 12 strips, we find X/j = 16. With an esti- 
mated 'socket' contribution of « 14.4 provided by surface 
evolution, we attribute a fraction of Xyk ~ ^-^ to statistical 
noise. This fraction decreases to « 1 if we use 30 strips, in 
which case we obtain a longitudinal resolution of 12°, com- 
parable to that achieved by Lanza et al. (2009). According 
to our test runs, we obtain reasonably stable results us- 
ing 12 strips. As the stability of the solutions decreases for 
larger strip numbers, while only slightly improves, we 
argue in favor of using 12 global strips in our modeling, to 
extract the largest possible amount of physically relevant 
results. 

The resolution used on the eclipsed surface band, is de- 
termined according to the following considerations. The ex- 
tent of the planetary disk on the center of the stellar disk is 
about 20° X 20° . All stellar surface elements simultaneously 
(un)covered by the planet's disk are equivalent with respect 
to our lightcurve modeling. Individual features can, thus, 
be located (or smeared out) all along the edge of the plane- 
tary disk to provide the same effect in the lightcurve. This 
edge stretches across 10° in longitude (only the 'forward' 
part) and 20° in latitude, which defines a fundamental limit 
for the resolution. Assuming a particular shape for the fea- 
tures, decreases the degree of ambiguity as was for example 
shown by Wolter et al. (2009). 

A meaningful structure in the transit profile should 
comprise at least 3 consecutive lightcurve bins correspond- 
ing to about 360 s or « 6° of planet movement across the 
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center of the stellar disk. The extent of individual strips 
should, therefore, not fall below this limit; it should even 
be larger. 

Combining these arguments with the results of our test 
runs, we decided to use 24 strips on the eclipsed section, so 
that a longitudinal resolution of 15° is achieved. With this 
choice, a single strip on the eclipsed band appears about 
the same size (face-on) as the planetary disk. Additionally, 
we note that this approximately corresponds to the resolu- 
tion used by Lanza et al. (2009) in their maximum entropy 
reconstructions . 



3. 6. Results of the modeling 

In our analysis, we achieve a longitudinal resolution of « 15° 
on the eclipsed section making up « 17 % of the stellar 
surface and 30° for the rest. 

In Figs. 3 and 4 we present the results of our modeling. 
Figure 3 shows the entire sub-sample of CoRoT data points 
used in our modeling as well as our lightcurve model in the 
upper panel. In the lower panel we show the model residuals 
(see Sect. 3 for the definition of the error). Obviously, the 
data are matched well; however, there are systematic offsets 
between the observation and the model. In particular, the 
model tends to overestimate the observations during the 
first half of the time span, whereas it underpredicts it in 
the second half. This eS'ect is related to surface evolution 
already detectable on time scales below the rotation period 
(Lanza et al. 2009) (also see Sect. 3.4). Within the tran- 
sits the residuals remain small compared to the rest of the 
lightcurve. This must be regarded a consequence of both 
the smaller bin size of 128 s used here and also the twofold 
better resolution of the model on the eclipsed section. Note 
that during the fit the transit residuals are 'overweighted' 
by a factor of ten to avoid them to be prevailed by the much 
larger global residuals. Although, the deviations can be as 
large as 10 a, the mean deviation of the global lightcurve 
from the model amounts to 620 x 10~^ not far from the 
theoretical limit of 560 x 10"^ (cf. Sect. 3.4). 

The lightcurve presented in Fig. 3 contains six tran- 
sits (labeled 'Tl — 6'). The associated transit lightcurves 
together with our models are shown in detail in Fig. 4. 
Each individual panel shows the same transit twice: The 
lower curve represents a transit reconstruction from the full 
data sample (phases 1.85 — 3.85), and the upper curve de- 
notes a reconstruction from only the first (Tl — 3) or second 
(T4 - 6) half of the sample data (shifted up by 0.03). The 
dotted lines show the transits as we would observe them 
without any activity on the eclipsed section of the surface, 
where we assume a brightness of 1 for the underlying pho- 
tosphere. 

The transit reconstructions obtained from half of the 
sample data reproduce the transit substructure very accu- 
rately. The resulting surface reconstructions are, however, 
unreliable where the surface is insufficiently covered (cf. 
Fig. 2, around longitudes of 180° and 320°). Interestingly, 
those reconstructions based on data from two rotation 
phases also recover most of the transit substructure and 
are by no means off the mark. When both rotation periods 
are used, typically increases by 10 — 20 %, a difference 
hardly visible in Fig. 4. As an exception, the fit quality of 
the third transit (T3), decreases dramatically, with ™" 
creasing by a factor of « 2.5. This is, however, mainly a 
consequence of the observed surface evolution shifting the 



continuum level. The overall stability of the fit quality indi- 
cates that lifetimes of surface features are of the order of a 
few stellar rotation, which is in agreement with the results 
of Lanza et al. (2009). 

In Fig. 5 (lower and middle panel) we present the strip 
brightness distribution pertaining to the lightcurve model 
shown in Figs. 3 and 4, i.e., a ID-reconstruction of the 
surface. We estimated mean and errors by recording the 
distribution of the parameter values obtained from 50 re- 
constructions with randomized starting points, and the re- 
spective distributions are indicated by the color gradients in 
Fig. 5. The error bars correspond to the associated standard 
deviations. They refiect the ability of the fitting algorithm 
to converge to a unique extremum, which is determined by 
both the characteristics of the algorithm and the structure 
of the fit statistics. Investigating the brightnesses distribu- 
tion of the non-eclipsed strips, we notice a slight degeneracy 
in some of the 12 non-eclipsed strips, i.e., a fraction of the 
brightness may be redistributed without considerable loss 
of fit quality. The averaging of the 50 reconstructions fiat- 
tens out such features, and, thus, acts like a regularization 
of the brightness distribution. No such effect is observed for 
the eclipsed strips. 

We compared our results to the reconstructions given 
by Lanza et al. (2009) (their Fig. 4) and find our longitude 
scale to be shifted by sa 70° in respect to the Lanza et al. 
scale. Our reconstructions show the same bright band at a 
longitude of ~ 260° (330° in our work) . Tentatively averag- 
ing over an appropriate 'time band' in their Fig. 4, we also 
find qualitative agreement for the remaining spot distribu- 
tion. 

Clearly, the flux fraction contributed by the eclipsed 
strips is smaller than that of the non-eclipsed strips, be- 
cause they are smaller by a factor of five. In the upper 
panel of Fig. 5 we show the lightcurve model contributions 
provided by the eclipsed and the non-eclipsed section with 
their sum making up the model for the CoRoT data which 
is also shown. The median flux level was subtracted from 
all curves to emphasize the modulation amplitude in favor 
of flux level. Obviously, the modulation amplitudes induced 
by the eclipsed and non-eclipsed section approximately bal- 



MJD [dl 



0.98 

0.97 

0.96 

0.95 

0.94 

0.93 

0.92 

0.91 

0.9 
' 10 


-10 




Fig. 3. Upper panel: CoRoT data for rotational phases 1.85 to 
3.85 (2 016 s binning for the global lightcurve and 128 s for the 
transits, red symbols) and our model lightcurve (dashed blue 
curve). Lower panel: Residuals of our model. 
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Fig. 4. Close-up on the six individual transits (128 s binning). 
Observations are drawn as red points (including la errors). The 
upper curve (black) in each panel shows the fit obtained from 
1 rotation (phases 1.85 to 2.85 for Tl - 3 and 2.85 to 3.85 for 
T4 — 6); the lower curve (blue) gives the fit obtained by modeling 
both rotations (as seen in Fig. 3). The dashed lines show the 
undisturbed transit profile for comparison. 



ance. This indicates that their influence on the stellar vari- 
ability is of the same order despite their large difference in 
size. 

A visualization of our best fit surface reconstruction is 
presented in Fig. 6, showing a Hammer projection of the 
associated distribution of surface brightnesses. The planet- 
defined low latitude band shows especially dark features e.g. 
in the range of 200° to 300° in longitude and is clearly visi- 
ble. Also the non-eclipsed sections of the star show signifi- 
cant variations. Note that our map only shows the average 
brightness of these regions; since the non-eclipsed regions 
are larger by area, they contribute more flux, however, the 
'missing' flux in these regions is likely also concentrated in 
spots. In the following section we address the issue of the 
flux contribution from the eclipsed and non-eclipsed sec- 
tions. 



3.6.1. Brightness distribution and spot coverage 

Without a very precise absolute flux calibration (as e.g. 
in Jeffers et al. 2006), lightcurve analyses can usually only 
investigate the inhomogeneous part of the entire spot cov- 
erage. This statement is, however, partially invalidated by 
a transiting planet because it breaks the symmetry of the 
problem: Spots being eclipsed by the planetary disk dis- 
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Fig. 5. Upper panel: The CoRoT data (red symbols) and the 
fiux modulation contributed in our model by the eclipsed (blue 
curve) and non-eclipsed section (dotted gray curve) individu- 
ally. All curves are median subtracted. Middle and lower panel: 
Brightness distribution of the strips located on the eclipsed and 
non-eclipsed section of the surface. The color gradient renders 
the distribution obtained from 50 reconstructions with random- 
ized starting points. 




Fig. 6. Surface map of CoRoT-2a showing the reconstructed 
brightness distribution. Spots located on the non-eclipsed sur- 
face are, due to their unknown latitude, blurred over the entire 
reconstruction strip resulting in the lower contrast compared to 
the eclipsed section. 



tort the transit proflles no matter whether they belong to 
a structure which appears symmetric on a global scale or 
not. 

As an example, assume that half the eclipsed section 
of CoRoT-2a, say longitudes — 180°, is spotted while the 
other half is covered by undisturbed photosphere. Clearly, 
the transits will be shallower when the planet eclipses the 
dark portion of the star, and they will be deeper when the 
bright section is eclipsed. Also the lightcurve will be dis- 
torted. Now further assume that a comparable section be- 
tween longitudes 180 — 360° is dark on the opposite hemi- 
sphere of the star outside the eclipsed band. In this case the 
global spot conflguration is perfectly symmetric with re- 
spect to longitude and the global lightcurve does not show 
any trace of activity. Yet, the transits will still be shallower 
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when the planet echpses the dark band, and they will still 
be deeper when the bright surface is eclipsed. 

If the spots were distributed symmetrically across the 
stellar surface, we would expect the stellar surface to be 
homogeneously bright. In Fig. 7 we show the brightness 
ratio of eclipsed and non-eclipsed section as a function of 
longitude. Since there are more eclipsed than non-eclipsed 
strips, we always compare strips covering the same longi- 
tude. Only in two cases the eclipsed section is brighter than 
its non-eclipsed counterpart, while in 22 cases it is not. 

The mean ratio is 0.94±0.01 so that the part of the star 
passingly covered by the planet is found to be 6 % darker 
than the rest of the surface. Note that the remaining (non- 
transited) surface is brighter on average, locally it may even 
be darker. 
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Fig. 7. Ratio of the surface brightness in the eclipsed {Beci. ) 
and non-eclipsed {Bnon-eci.) section. 



4. Discussion and conclusions 

We present a surface reconstruction for the planet host star 
CoRoT-2a. Our modeling is based on a CoRoT data in- 
terval covering two full stellar rotations, and it treats the 
entire lightcurve - including the planetary transits - in a 
consistent way. 

We show that a consistent modeling of the lightcurve is 
possible using a static model, i.e., not including any spot 
evolution. Although surface evolution on scales of the stel- 
lar rotation period is seen in both the reconstruction of the 
global lightcurve (as already reported by Lanza et al. 2009) 
and the transit lightcurves, this effect is small in the context 
of our analysis. The static model provides reasonable fits to 
six consecutive transit lightcurves. The associated surface 
configuration changes little during this period, and, there- 
fore, the surface evolution must be relatively slow compared 
to the time scale of « 9 d under consideration. This time 
scale is also valid for the lifetimes of spots on the eclipsed 
surface section. 

Our results indicate that the planet-eclipsed band on 
the stellar surface is - on average - about 6 % darker then 
the remaining part of the surface. Lanza et al. (2009) note 
that the strength of differential rotation derived from their 
lightcurve fits seems very low in comparison to expected 
values derived from measurements in other systems (Barnes 
et al. 2005). They speculate that this may indicate a spot 
distribution limited to a narrow latitude band. If this should 
be true, the latitude band is possibly located at low lati- 
tudes, i.e., within ±30° around the equator as observed on 
the Sun. In this case it covers the eclipsed section where 
we find a darker surface, i.e., higher spot coverage. We cau- 
tion that this result may also be influenced by the adopted 
planetary parameters (mainly the size), which are hard to 
determine accurately (Czesla et al. 2009). 



We checked whether the effect of gravitational dark- 
ening could significantly contribute to a darker surface in 
the vicinity of the equator. For the stellar parameters of 
CoRoT-2a, we find that the (effective) gravitational accel- 
erations at the poles and at the equator are equal to within 
0.07 %, so that gravitational darkening does not provide a 
significant contribution to the brightness gradient found in 
our modeling; this result is nearly independent of the as- 
sumed coefficient, /3i {T^jf « g^^ with the effective surface 
temperature Tk// and the surface gravity g), which is of the 
order of 0.3 - 0.4 for CoRoT-2a (Claret 2000). 

The 'narrow-band hypothesis', i.e., a higher spot cov- 
erage in the planet-eclipsed section compared to the non- 
eclipsed surface, also provides a natural explanation for the 
fact that both the eclipsed and non-eclipsed surface regions 
account for about the same amplitude of variation in the 
lightcurve. Using the Sun as an analogy again, we would 
qualitatively expect the same structure, as seen under the 
planet path, on the opposite hemisphere as well; two 'active 
belts' that are symmetric with respect to the equator. The 
non-eclipsed activity belt, which would be only observable 
in the global lightcurve, would then be primarily responsi- 
ble for the variability of the lightcurve contributed by the 
non-eclipsed surface section. 

We conclude that our results support a surface model 
consisting of active regions north and south of the equa- 
tor, possibly even bands of spots at low latitudes analo- 
gous to the Sun. Further investigations of this system us- 
ing more sophisticated models (first of all surface evolution) 
and using the entire observation interval of approximately 
140 days have the potential to reveal more information on 
the constantly changing surface distribution of spots on 
CoRoT-2a. 
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